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Abstract 

Using a sample of 25.9 million tp(2S) decays acquired with the CLEO-c detector at the CESR 
e + e~ collider, we report branching fractions for the decays XcJ — > ppn , pprj and ppu:, with 
J = 0, 1, 2. Our results for B(x c j — > ppn ) and B(x c j — > PPV) are consistent with, but more precise 
than, previous measurements. Furthermore, we include the first measurement of B(xcj ppu). 
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Recent theoretical results jlH3(] have highlighted the value of studying \l> — > ppM hadronic 
decay processes, where \l/ represents any cc bound state and M is a light meson. The 
application of these models allows measured ^ — > ppM partial widths to be used to estimate 
the production cross sections for a(pp — > \[>M), circumventing the calculation of some of the 
complicated underlying QCD processes. Calculations of this sort are interesting, for example, 
in the context of the future PANDA experiment [IJ which will exploit associated charmonium 
production in pp annihilation (pp — > \&M) in its search for exotic charmonia. Since the 
values of — > ppM) serve as key inputs for these calculations, the same authors have also 
developed techniques for calculating T(^/ — > ppM) j^], which can be tested with experimental 
data. In their meson emission model, they assume the sequential decay — Y pp — >■ ppM, and 
by applying techniques developed in (ij-yj, they estimate T(^f — >■ ppM) using the measured 
^ — > pp widths and well-known ppM coupling constants. If this sequential decay mechanism 
is in fact the dominant means by which ^ — > ppM decays proceed, then the branching 
fractions to ppM final states would provide a means of extracting other meson-nucleon 
coupling constants 

This article describes measurements of the branching fractions for Xcj decays to three 
final states, ppn°, ppr], and ppu using a sample of 25.9 million ip(2S) decays produced in 
e + e~ collisions at the Cornell Electron Storage Ring (CESR). The first observations of Xcj 
decays to ppn° and ppr\ were made by CLEO using a substantially smaller dataset of 3 million 
ip{2S) decays j6[. The ip{2S) produces a copious number of Xcj mesons via its radiative El 
transitions, ip(2S) — » 7Xcj, with branching fractions of approximately 9% for each of J = 0, 
1, and 2. We fully reconstruct decays of these secondary charmonia into pp'j'j and pp7i + ii^ii 
final states using the CLEO-c apparatus. 

The nearly hermetic CLEO-c 0] detector covers 93% of the solid angle. It features 
a 1 T superconducting solenoid housing drift chambers for tracking and charged particle 
identification and a ring imaging Cerenkov (RICH) system to further differentiate between 
charged particle species. Also within the solenoid volume is an electromagnetic calorimeter 
composed of 7784 CsI(Tl) crystals. The photon energy resolution in the calorimeter is 
2.2% at 1 GeV and 5% at 100 MeV and the momentum resolution achieved using the drift 
chambers is typically 0.6% at 1 GeV/c. 

In this analysis, we select events with either two or four charged tracks and at least three 
photons. Candidate tracks are required to have momentum p > 18.4 MeV/c and originate 
within a 10 cm long, 2 cm radius cylindrical volume centered around the e + e~ interac- 
tion point. The n^ 1 candidate tracks are required to have specific ionization measurements 
(dE/dx) consistent with those expected for charged pions within 3 standard deviations. The 
proton and anti-proton candidate tracks are required to have dE/dx measurements within 4 
standard deviations of the expected ionization losses for protons and anti-protons. If RICH 
information is available for the event, it is used in conjunction with dE/dx information to 
form joint likelihoods based on the hypothesis that the track is a proton, pion or kaon. 
Candidate protons and anti-protons are then required to be more proton-like than pion-like 
or kaon-like. Candidate photons are identified via the electromagnetic showers produced 
when incident on the calorimeter. They must be associated with showers depositing more 
than 30 MeV and have good separation from charged tracks. Neutral pions and rj mesons 
are then reconstructed in their 77 decay modes. The invariant mass of each photon pair is 
calculated, and the pair is accepted as a neutral meson decay candidate when the invariant 
mass is within 3 standard deviations of the corresponding meson's rest mass. The detected 
four-momenta of all final state particles are improved via a series of kinematic fits. In the 
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an d pp7i + 7i^7T° modes, the four-momenta of the photon pairs in 7r° and 77 candidates 
are constrained to the respective nominal rest masses taken from the Particle Data Group 
(PDG) report jsj. All charged tracks are constrained to originate from a single event vertex 
that is within the beam spot, which is measured on a run- by-run basis by tracking the typical 
location of the event vertex. Finally, the El photon, p, p and 7r°, 77 or 7r + 7r~7r° four-momenta 
are kinematically constrained to the initial state ip(2S) four momentum. Events are then 
selected according to the x 2 /d.o.f. of this four-constraint fit. For pp r j r ) final states, events 
are required to satisfy x 2 /d.o.f. < 5 and for pp7i + n^7T , x 2 /d.o.f. < 10. In instances where 
there is more than one possible combination of final state particles, e.g., more than one pair 
of photons satisfy the 7r° selection criteria, the four constraint kinematic fit is performed for 
all possible permutations, and the combination with the lowest xVd.o.f. is selected. This 
ambiguity in the final state particles occurs in around 10% of events. 

Further selection criteria to suppress backgrounds were investigated using a Monte Carlo 
(MC) sample of 1 x 10 8 ip(2S) decays generated using known partial widths from the PDG [s| 
and the models described in Ref. [9J for any unknown branching fractions. The dominant 
background in the ppn° mode is from ip(2S) — > 7i°7i°J/ifj, J/ifi — >■ pp events in which one 
photon from the two 7r° decays is soft enough not to skew the four-momentum to the extent 
that the event would fail the xVd-o.f. requirement. To suppress these, we reject events 
when the invariant mass of the pp system is close to the J /if} mass, that is, when 3.07 < 
M(pp) < 3.14 GeV/c 2 . The total remaining background does not peak in M(pp7c°) near 
any of the Xcj masses and accounts for 6% of the data passing our selection criteria. There 
is another small background in the ppr][^\ mode which peaks in M(ppr]) at the Xc2 mass 
from the process: Xc2 7^/Vs J/i> — > 1PP- The rate of this background is reduced by 
the initial requirement that M{ r y r y) is close to the 77 rest mass, and it is further suppressed 
by rejecting events when the invariant mass of the jpp system is close to the J/ip mass, 
i.e., when 3.07 < M(jpp) < 3.12 GeV/c 2 . The remaining background accounts for 30% 
of selected events but has no structure in M(ppr]). Additional requirements on kinematic 
variables are ineffective in suppressing background in the pp7r + 7r _ 7r° final state since this 
is dominated by Xcj ppit + n~ it decays in which the three pions do not result from an 
77 or bj. These events conserve four-momentum and cannot be distinguished from resonant 
(r],u) 7r + 7r~7T° decays on an event-by-event basis. Instead, a term is included to account 
for this background in a fit to the rj- and w-components of the decay. 

The branching fractions for Xcj — > ppM are calculated according to 

B (Xcj -)> PPM) = Nm . (1) 
eMN^2S)B 7 jB Y 

A^,(2S) is the number of if>(2S) present in the data [lOj. The signal efficiency of the combined 
CLEO-c apparatus, reconstruction and event selection algorithms, €m, is evaluated via anal- 
ysis of MC samples. The branching fractions for ip(2S) •yXcj, B 7 j, are those measured by 
CLEO Values for By, which represents the branchingfractions for M — > Y, where Y 
represents either 77 or 7r + 7r~7r°, are taken from the PDG [8|. The signal yield, Nm, is ob- 
tained via unbinned maximum likelihood fits to the data using slightly different techniques 
for the two distinct final states ppyy and pp7i + 7i^7i . 

In the first case, the yield is extracted via separate one-dimensional unbinned extended 
maximum likelihood fits to either the M(pp7c°) or M{jpprj) spectrum. A linear background 
term is included in both fits to account for the small, flat background surviving our selec- 
tion criteria. The signal shapes are modeled by Breit-Wigner distributions convolved with 
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FIG. 1: Candidate Xcj mass spectrum for Xcj — > ppn (a) and Xcj —> Ppflill] (b). Points with 
error bars are data, the solid lines show the fitted functions, and the dashed lines represent the 
linear background components of the fits. 

Gaussian resolution functions. The masses and widths of the Breit-Wigner distributions are 
fixed at the PDG values jsj for the Xcj- The resolutions are fixed at values extracted from 
MC simulations. The results of these fits are shown in Fig. [TJ 

The large non-77, non-o; background in the pp7i + 7i^7i° channel led us to choose a signal 
extraction method consisting of a two-dimensional unbinned extended maximum likelihood 
fit in M(pp7r + 7T~7T°) and M(7r + 7r~7r°) to simultaneously extract all six desired yields. Fitting 
in both variables provides sensitivity to the non-resonant background shape over a wide 
range of M(7r + 7r~7r°), which allows for a precise determination of the contribution in the r) 
and u signal regions. The ppu signal shapes are modeled as the product of Breit-Wigner 
distributions centered at the Xcj masses convolved with Gaussian resolution functions in 
M(ppiT + 7i^7T°) and a Breit-Wigner function centered at the w-mass and convolved with a 
Gaussian in M(tt + 7t~7t°). Similarly, the ppr] signals are modeled as the product of Breit- 
Wigner distributions convolved with Gaussian resolution functions in M(p/J7r + 7r~7r°) and, 
since the 77 is sufficiently narrow, a Gaussian in M(-k + 7t~tt ). The non-resonant background is 
represented by a reversed ARGUS function [l2| in M(7r + 7r~7r°) multiplied by the convolution 
of Breit-Wigner functions and Gaussian resolution functions in M(pp7i + 7i^ii ). The ARGUS 
function threshold is fixed at the three pion mass threshold. The non-peaking background 
can be well-described by a function linear in both M (ppir + -k~ -k°) and M (7r + 7r~7r°). The 
masses and widths of the Breit-Wigner distributions are fixed to the PDG values |8j. The 
Gaussian resolutions are allowed to float, although this parameter was constrained to be the 
same for each of the Xcj's. The probability density function (PDF) that provides the best 
fit to data is shown in Fig. [5J and projections of the data and fit onto the M(pp7r + 7r~7r°) 
and M(tt + 7t~7t°) axes are shown in Figs. [3]and|H 

The branching fractions are assigned systematic errors due to uncertainties in charged 
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3: Results of the two-dimensional fit to M(7r + 7r _ 7r°) and M(ppir + ir~iT ) for the 
ppTr + ir~7t° channel projected onto the Xcj candidate mass axis for three M(7r + 7r _ 7r°) re- 
gions: (a) 0.541-0.553 GeV/c 2 , (b) 0.770-0.794 GeV/c 2 , (c) 0.560-0.730 GeV/c 2 . Points with error 
bars are data, the solid lines are projections of the total fitted function, the dashed lines are 
the summed background components of the fit, and the dotted lines show the planar background 
components. 
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FIG. 4: Results of the two-dimensional fit to M(7r + 7r~7r°) and M(ppir + ir~iT ) for the 
Xcj ppTT + ir~'K channel projected over the entire range of M (ppTr + ir~ ir°) onto the M(7r + 7r _ 7r ) 
axis. Points with error bars are data, the solid line is a projection of the total fitted PDF, the 
dashed line is the summed background components of the fit, and the dotted line shows the planar 
background component. 



particle tracking efficiency (1% per track), uncertainties in photon finding efficiency (2% 
per photon), the uncertainty in the number of ip(2S) (2%) [To| and an error propagated 
from the uncertainty in the ip(2S) — > •"/Xcj branching fractions (5%-7%) The x 2 /d-o.f. 
requirement introduces a small uncertainty (0%-3%) in the selection efficiency which is 
assessed by varying the requirement and repeating the analysis. The error introduced from 
the choice of fitting technique is estimated by varying the fitted PDF and observing the 
change in efficiency-corrected yield. Specifically, for the pp'yy final states, the Xcj masses 
are varied by one standard deviation from the nominal PDG values [§[ and the systematic 
uncertainty assigned as the average variation in observed yield between the two extremes 
of Xcj mass. A similar variation is repeated for the Xcj widths. The linear background 
term is replaced with a second order polynomial and the error assigned as the change in 
observed yield. This uncertainty associated with the parametrization of the background 
shape is the dominant contribution to the fitting systematic and is largest in the ppr] final 
state (l%-8%). The experimental resolution extracted from MC events is varied by ±25% 
and the error assigned as the average yield for the two extremes. The errors from each of 
these individual variations are added in quadrature to obtain the total error associated with 
the fitting technique. 

A similar set of variations is applied to the two-dimensional fits to the pp7r + 7r~7r° final 
state. The Xcj masses and widths are varied by one standard deviation from the PDG 
values jsj. This is repeated for the mass and width of the u and the mass of the rj. In 
addition, the ARGUS function is replaced with a third order polynomial and the change in 
observed yield assigned as the uncertainty. Finally, the range of M(tt + -k~7t ) fitted is varied 
from [0.41,1.0] to [0.53,0.85] and the uncertainty taken as the change in observed yield. 
Again, the individual errors are summed in quadrature to obtain the total error associated 
with the fitting technique (l%-5%). 

The efficiencies are extracted from MC simulations in which the Xcj ppM decays pop- 
ulate phase space uniformly; 6m is obtained by fitting either the ID M{ppM) distribution or 
the 2D M(ppM), M(pp7i + 7i^7i°) distribution. As a result, €m is an average of the efficiency 
over the ppM Dalitz plot. We have investigated both how the data populate the Dalitz plot 
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FIG. 5: Dalitz plots for the ppir final state in decays of (a) XcO an d (b) Xci- The contours indicate 
the density distributions predicted by the meson emission model [jj. 

and how the efficiency varies across the Dalitz plot. Initially, Xcj candidates are selected in 
data via a requirement on M(ppM), and the surviving events are binned in terms of M 2 (pM) 
and M 2 (pM). The unbinned Dalitz plot is shown in Fig. |5]for Xc0,i — > pp7c°. It is clear that 
the data do not populate phase space uniformly. In particular, a broad structure close to pp 
threshold is evident in Xco — > PP^°- The same selection criteria are then applied to the signal 
MC, and again the MC sample is binned in terms of M 2 (pM) and M 2 (pM). By dividing 
out the number of generated MC events in each bin, we obtain the efficiency as a function of 
M 2 (pM) and M 2 (pM). To obtain the efficiency as a function of M 2 (pM) and M 2 (pM) for 
the Xcj — > Ppu channels, a sideband subtraction in M(n + n^n ) is also required to suppress 
the non-resonant background (see Fig. HJ). We find that the efficiencies are smoothly varying 
in M 2 (pM) and M 2 (pM) in all channels with the exception of the region of M(pp) close to 
the J/ijj mass in the Xcj — > ppvr channels. In this region the efficiency is considerably lower 
as a result of the requirement on M{pp) needed to suppress backgrounds. 

Variations in detection efficiency across the Dalitz plot together with resonant structures 
in the data could potentially lead to a large systematic uncertainty in the value of e*/. We 
quantify this uncertainty by calculating the efficiency corrected yield in two ways. First, 
the Dalitz plot for the data are integrated to obtain the yield and this is corrected using 
the efficiency averaged over M 2 (pM) and M 2 (pM), a procedure that gives an efficiency- 
corrected yield close to that used in the nominal analysis. Next, the data are corrected 
for efficiency as a function of M 2 (pM) and M 2 (pM). The systematic uncertainty is then 
assigned as the difference in efficiency corrected yields obtained using these two methods 
and is less than or equal to 5% in all channels excluding Xco ppf][7i + ii^'n- ]. In that case 
the uncertainty is slightly higher at 9%. 

The final product branching fractions, B^j x B(x c j ppM) , are given in Table |T] where 
the errors shown are statistical then systematic. The final branching fractions for the pro- 
cesses Xcj ppir , ppp and ppu are given in Table UTI In this case the errors are statistical, 
systematic due to detector and analysis uncertainties, and a separate systematic error due 
to the uncertainty in the ip(2S) — > 7% c j branching fractions. A weighted average of the two 
separate B(x c j ppp) measurements is made taking into account correlated systematic er- 
rors. The results are in good agreement with the previously published CLEO data [6j with, 
as expected, a factor of around 3 smaller statistical errors. 

The meson emission model predictions for the two branching fractions 
B{Xco -> PPTT°)thcory = 2.5 x 10~ 4 and B(xd -> PFr°)theor y = 0.2 x 1(T 4 jl] are well be- 



low our observed branching fractions, by factors of about 3 and 10 respectively. This 
suggests that meson emission, as described by this model, is not the dominant decay 
mechanism. This can be further demonstrated by comparing the theoretical Dalitz plot 
event densities calculated in reference jsj with our data; this comparison is shown in Fig. 
The meson emission model predicts strength in regions of low pir and p7r° invariant mass, 
whereas the data show a clear enhancement at low pp invariant mass. 

In summary, we have presented new measurements of the branching fractions 
B(Xcj — » PpK ), B(xcj — > PPv) an d &(Xcj —> Ppw). We find good agreement with the avail- 
able previous experimental results for Xcj — > ppn° and Xcj — > PPV, and the large CLEO 
data set allows us for the first time to disentangle the XcJ — > PP^ strength from the large 
non-resonant background. Finally, we make a comparison with the meson emission model 
calculations of Barnes et al. and find that the sequential emission process described by the 
authors does not describe our data. 
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TABLE I: The product branching fractions B = B^j x B(\ c j — > ppM). Uncertainties are statistical, then systematic. 6m is the signal 
efficiency (described in the text). Yield is the number of signal events evaluated via an unbinned maximum likelihood fit (described in the 
text). 



XcO Xcl Xc2 





6M(%) 


Yield 


£(10" 5 ) 


e M (%) 


Yield 


£(10- 5 ) 


e M (%) 


Yield 


£(io- 5 ) 


ppir 


33.4 


610.9 


7.15 ±0.34 ±0.47 


36.1 


146.9 


1.59 ±0.15 ±0.12 


35.2 


406.6 


4.51 ±0.24 ±0.33 




30.6 


99.4 


3.18 ±0.42 ±0.23 


33.1 


49.9 


1.48 ±0.25 ±0.17 


33.1 


58.1 


1.72 ±0.27 ±0.16 


pp7][ir + 1T~ IT ] 


21.3 


51.9 


4.15 ±0.61 ±0.47 


23.4 


18.1 


1.33 ±0.33 ±0.12 


22.5 


19.9 


1.51 ±0.37 ±0.12 


PpLO 


22.2 


263.0 


3.43 ±0.35 ±0.26 


23.9 


113.7 


1.42 ±0.20 ±0.13 


23.4 


185.5 


1.64 ±0.22 ±0.13 



TABLE II: Final Xcj — > ppM branching fractions. Uncertainties are statistical, then systematic and then a separate systematic error due 
to the uncertainty in the ip(2S) — > JXcJ branching fractions. 





J = 


B (xcj ^PPM) (10^) 
J = 1 


J = 2 


PpTT° 

pp^fmean] 
ppto 


7.76 ±0.37 ±0.51 ±0.39 
3.73 ±0.38 ±0.28 ±0.19 
5.57 ±0.48 ±0.42 ±0.28 


1.75 ±0.16 ±0.13 ±0.11 
1.56 ±0.22 ±0.14 ±0.10 
2.28 ±0.28 ±0.16 ±0.14 


4.83 ±0.25 ±0.35 ±0.31 
1.76 ±0.23 ±0.14 ±0.11 
3.68 ±0.35 ±0.26 ±0.24 



